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The main part of this talk is a review and summary of how QCD is used in two main 

areas of nuclear physics, namely in determining the quark flavor and spin content of the 

^ ' proton and in ultrarelativistic heavy ion collisions. Brief comments are made concerning 

^ ■ effective theories in hadron physics and on the separation of various twists in using the 

Q ! operator product expansion to analyze hard processes. 
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0\ • Quantum Chromodynamics (QCD) is a fundamental theory in the sense that one ex- 

r-| ■ pects QCD to exist not just as a perturbative expansion but also in its strong coupling 
5^. regime. There is good evidence from lattice QCD calculations that this is indeed the case. 
Qh! QCD stands in contrast to QED or the Electroweak Theory which have good perturbative 
(-H I expansions but which are not expected to exist in the strong coupling regime. In partic- 
ular it is widely believed that the Electroweak Theory must be unified with additional 
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1. Introduction 



interactions at a scale below its Landau ghost. 

The likelihood that only asymptotically free theories make sense puts severe restrictions 
5^ I on the use of "effective" field theories, field theories that are to be used in a limited range of 
scales. For example, it probably does not make much sense to try and describe low energy 
pion-nucleon interactions in terms of a pseudoscalar meson-nucleon interaction Lagrangian 
since the coupling necessary in such a description is so strong that the theory is internally 
inconsistent even at low energy scales. Of course the Born terms (tree graphs) of such 
a theory still make sense since they represent an analytic structure which is completely 
determined by the lowest lying states of the pion-nucleon system. 

There is a special property of pion interactions which makes it possible to use particular 
effective theories as a description of low energy pion-nucleon interactions. That special 
property is chiral symmetry. One of the properties of the pion, as the Nambu-Goldstone 
boson of chiral symmetry breaking, is the fact that its couplings become weak at low 
momentum. Thus if one uses an effective chiral theory to describe low energy pion- 
nucleon interactions, and if a cutoff is introduced so that the couplings never grows large, 
then one can use that effective theory, even in a nonperturbative manner, to make precise 
calculations [|l]-|^. Such theories tend to have quite a few parameters since the cutoff is 
not very high and all the dynamics above the cutoff must be put into parameters of the 
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effective tlieory. Tliere lias been mucli interest lUg] and progress in tlie last few years 
using effective theories to describe light nuclei in terms of pion-pion and pion-nucleon 
interactions. It is perhaps a bit too early to decide whether this description is superior 
to the traditional potential approach, but one can hope to get better insight as to where 
it is more useful to use quarks and gluons to describe hadronic physics and where it is 
better to use effective degrees of freedom. 

As far as is known Q CD- type theories are the only four-dimensional field theories 
"known" to exist as fundamental theories. The goal in strong interaction physics is to 
understand how QCD works. Few physicists doubt the validity of QCD as the theory of 
the strong interactions. However, QCD is a rich and sophisticated theory about which 
there is still much to understand and appreciate. 

I shall cover three topics in this talk. Firstly, 1 shall review progress in understanding 
the quark and gluon structure of the proton wavefunction. This has become one of the 
most active fields in nuclear physics in which the Jefferson Laboratory has now begun 
to make important new contributions. This part of my talk will have some overlap with 
and is complementary to the talks of L. Cardman|0 and A. MagnonJ^]. The next topic 
to be covered is that of heavy ion collisions and small-x physics. The object here is 
to describe the relationships between the physics being pursued in high energy small-x 
physics and that which will become available at heavy ion colliders at Brookhaven and at 
CERN. The final topic is perhaps more of a remark, but it is a remark which is important 
for medium energy physics. The point is that the separation between perturbative and 
nonperturbative QCD is not always so easy to define. Nevertheless, one must often make 
such a separation in medium energy phenomenology and care must be taken to insure that 
that the separation, even if "scheme" dependent, is done consistently. A good example is 
a recent analysis by Kataev0 and his collaborators of the description of the F3 structure 
function in QCD. 

2. Spin and Strangeness in the proton 

2.1. Spin and the constituent quark model 

In the constituent quark model the proton is made of two up-quarks and a down-quark. 
In the nonrelativistic version the three quarks are in zero orbital angular momentum 
states and so the spin of the proton is given by the sum of the spins of the quarks. With 
the usual SU(6) wavef unctions |llO[| , and with Ag labeling the z-component of spin of the 



quark q, one has 

AM = 4/3,Arf= -1/3 (1) 

leading to 

Ga = Au-Ad = 5/3 (2) 

and a total spin of the proton, AE, given by 

AS = An + Arf = 1. (3) 

The nonrelativistic quark model gives a reasonable picture of the proton but the value of 
Ga is clearly somewhat high. 



Relativistic quark models give better agreement with experiment for Ga- Typically ||11|| 
one has 

AM^l,Arf^-l/4 (4) 

giving 

Ga = ^u-M^ 5/4 (5) 

and 

AS = An + Arf ^ 3/4. (6) 

In neither the nonrelativistic nor in the relativistic versions of the model is there any room 
for strange quarks in the proton. In addition to Ga the quark model gives a very good 
account of baryon magnetic moments and it has been the basis on which spectroscopy 
of mesons and baryons has been discussed for some time[0. Much of what follows will 



be concerned with whether the quark model is able to give a reasonable account of the 
protons total spin. 

2.2. How to measure spin 

The axial vector current of flavor-f quarks is given by 

jL = ^/7M75g/- (7) 

For free quarks 

ips\J5f^\ps) = 2ms f, _^ 2\p^ (8) 

where s^ is the fermion spin four-vector corresponding to a quark with spin orientation 
s in its rest system, and where A is the quark helicity. In a frame where quarks have a 
large longitudinal momentum it is convenient to quantize spin along the direction of the 
large momentum. 

In the quark-parton picture of the proton it is natural to suppose that 

{Ps\jl^\Ps) = 2Mps,/\q} (9) 

where \Ps > is a proton state, and where Agj is the fraction of the proton's spin carried 
by quarks of flavor f . 

In spin-dependent deep inelastic scattering on a proton one can measure a particular 
combinarning: Citation 'lof ' on page 3 undefined on input line 179. ! Missintion of axial 
vector current matrix elements, that given by 



2Ms^ 



f^ dx gfix, Q') = li: e}{Ps\ji,\Ps). (10) 



Scattering on neutrons gives an independent combination of axial vector currents. Defin- 
ing 

/'dx,r^Hx,Q^)=rrHQ^) (11) 

Jo 

and using (9), one finds 

T^ = -(-Au + -Ad+-As) (12) 

and 

T^ = -(-Ad + -Au + -As). (13) 

^ 2^9 9 9 ^ ^ ^ 

An additional relation, involving Am, Ad and As comes from semi- leptonic hyperon 
decays|l2|. After using SU(3) flavor symmetry one finds 



3F-D = Au + Ad- 2As. (14) 

Equations (12) - (14) apparently (but see below) give enough information to determine 
Aqj from rf,rf,F and D. 

2.3. Bare versus constituent quarks 

In the constituent quark model the proton is described in terms of quarks which are 
not really point-like and which are not the quark degrees of freedom which appear in 
the QCD Lagrangian. The quark fields in (7) and (10) are the fundamental (bare) fields 
of the QCD Lagrangian. A basic assumption of the constituent quark model is that for 
static matrix elements one may replace the fundamental quark fields by the constituent 
(effective) quark fields of the quark model. This gives a method of calculating forward, or 
near forward, matrix elements of local bare quark currents in terms of wavefunctions of 
the constituent quark model. This is the basis on which one can see if there is agreement 
between the quark model and results obtained from deep inelastic scattering which do 
not directly measure constituent quarks but which, through relations like (10), do lead to 
static matrix elements of local currents. 

2.4. A final subtlety 

While (9) is a plausible assumption in fact it is not quite right. The correct relation is 



{Ps\ji^\Ps) = 2Ms,{Aqj - ^AG) (15) 

with AG the amount of the proton's spin carried by gluons. The argument for (15) is 
subtle and still not without controversy. I think the simplest argument for its validity 
comes from considering an example in quantum electrodynamics with massless electrons. 
At order a a virtual photon can break up into an electron-positron pair. Because helicity 
is conserved in electromagnetic interactions the helicities of the electron and positron 



will be of opposite sign. This if one takes the matrix element of js^ = ip'jfj.'j^ip between 
virtual photon states, at order a, one would expect to get zero since js^, from (8), is 
supposed to measure electron (and positron) helicities. (The situation is illustrated in 
Fig.l.) However, the matrix element turns out to be — ^, for a transversely polarized 
photon having helicity 1, despite the fact that any given electron-positron state gives a 
zero value. This is the axial anomaly and — ^ is the value of the axial anomaly. In QCD 
exactly the same phenomenon occurs. The — ^ in (15) is the anomaly value while AG 
gives the net helicity of gluons in the proton. The — ^ comes from effects which cannot 
be ascribed to the quark wavefunction of the proton and is one of the most profound 
elements in QCD. 




Figure 1. The axial vector current couples to electrons and positrons having opposite 
helicities 



2.5. Results from spin-dependent deep inelastic scattering 

There are several groups that have done extensive global fits, using a second order 
renormalization group formalism, to spin-dependent deep inelastic scattering |T7|-|19[ . The 
experimental situation is discussed in some detail in the talk of A. Magnon]^. Here, I 
shall briefly recount the results found by Altarelli et al who have carried out a series of fits 
to all the available data. The effect of polarized gluons is quite important in the various 
fits. The preferred fit has 



An = 0.86, Ad 
AG 



-0.37, As = -0.05 
1.4 ±0.9 



(16) 



giving 



AS = 0.44 ±0.09. 



(17) 



While these numbers are not in perfect agreement with the relativistic quark model the 
situation is not so bad either. In particular, the reasonably small value for As is much 



more comfortable for the quark model that were early values. However, because of the 
large errors on AG it is perhaps premature to take the final results as conclusive. The 
situation would be considerably improved by a direct measurement of AG, as planned at 
COMPASS and RHIC and which could also be done well at HERA. A few reliable points 
for AG(a;, Q"^) would help considerably to fix the first moment, AG, of AG(a;, Q^). 

2.6. Strangeness as viewed by vector currents 

In the previous sections we have seen how spin-dependent deep inelastic scattering 
gives information on the polarized quark and gluon distributions of the proton. New 
experiments at Bates and the Jefferson Laboratory give complimentary information on 
the number of strange quarks in the proton. By measuring very accurately parity violating 
elastic electron-proton scattering information on the contribution of strange quarks to the 
proton's electromagnetic form factors can be obtained. This is discussed in some detail 
in the talk of L. Cardman0 so I will simply state the results here and then comment on 
their significance and how they fit in with the general program of determining the quark 
and gluon content of the proton. If G|;/^,j-)(Q^) stands for the strange quark contribution 
to the electric (magnetic) form factor of the proton at a momentum transfer Q^ then the 



new result from the SAMPLE experiment at Bates ispO| 

Gl^iQ^ = O.lGeV^) = 0.23 ± 0.37 ± 0.15 ± 0.19 (18) 

where the first error is statistical, the second is a systematic error and the third is due to 
the axial form factor G^. The HAPPEX experiment at the Jefferson Laboratory [^ finds 



G% + 0.39 GIj = 0.023 ± 0.034 ± 0.022 ± 0.026 (19) 

at Q^ = 0A7GeV'^ where the first error is statistical, the second systematic and the third 
comes from uncertainties in the electric form factor of the neutron. To set the scale for 
these numbers we note that in the constituent quark model 

G^(0)=2.47,Gl,(0) = 0-32 (20) 

and 

G^(0)=4/3,G^(0) = -l/3 (21) 

give the up and down contributions to the proton's magnetic moment and charge, respec- 
tively. 

The HAPPEX result suggests that strange quarks give no more than a few percent of 
the up plus down contribution to Ge + 0.39Gm- The SAMPLE result is also consistent 
with a small strange quark contribution, but the limit is perhaps not too stringent. 

2.7. Summary on quark content of the proton 

Perhaps the main issue in determining the quark and gluon content of the proton is 
the issue of how many and "what kind" of strange quarks are to be found in the proton's 
wavefunction. At first sight this would seem to be a simple problem whose answer, at 
least roughly, has been known for some time. After all, spin-independent deep inelastic 



lepton-proton experiments have determined that x{s{x, Q'^) + s{x, Q"^)) is sizeable. Indeed, 
the strange quark sea is about one-half that of the non-strange sea. 

However, a moment's reflection is enough to realize that while it is interesting to know 
the number, and the Bjorken x distribution, of strange quarks in the proton it is even 
more important to know if those strange quarks are just short time fluctuations which 
play no role in the dynamics of the proton or if the strange quarks live long enough to be 
essential in determining the proton's mass and wavefunction. For example, the fluctuation 
of a gluon of the proton into an ss pair, in which the relative transverse momentum of 
the s and s is greater than a few GeV or so, contributes to the spin-independent strange 
sea distribution, but such pairs are too compact and short-lived to interact with the rest 
of the proton and so are uninteresting for static properties of the proton. In particular 
such short-lived fluctuations should not contribute to As, because the s and s helicities 
will cancel out in As, and they should not contribute to G\,j or G|;, again because the s 
and s will cancel. 

We can now begin to appreciate the new information contained in spin-dependent deep 
inelastic scattering and in strangeness as measured in parity violating elastic electron pro- 
ton scattering. As(Q^) measures the sum of the s and s helicity fractions in the proton 
independently of the longitudinal momentum fraction of the s and s and including all 
transverse momenta up to Q. (One expects no Q^-dependence of As so long as Q^ is grater 
than a few GeV"^.) In order that the s and s not have cancelling helicities the transverse 
momentum of the s and s must be small enough that they have helicity nonconserv- 
ing (nonperturbative) interactions in the proton or that the current quark mass, about 
ISOMeV, not be negligible. For example, an interaction correlating the s with a spec- 
tator u-quark, perhaps giving a virtual K^ meson, could certainly polarize the strange 
sea. Strange quarks measured by vector currents, as in the parity violating experiments, 
indicate a difference between the transverse momentum (or transverse coordinate) distri- 
butions of the s and s. This, again, requires that the s and/or the s interact with the 
remnants of the proton. 

Thus both As and G^ give a measure of strange quarks in the proton which have some 
interaction with the rest of the proton and hence are an integral part of the proton's wave- 
function. (There is, however, a contribution to As coming only from mass effects where 
interactions in the proton are not necessary.) Thus sizeable values of As and G"* would 
indicate that strange quarks play an essential role in the proton's wavefunction. If both 
As, after subtracting the anomaly, and G^ turn out to be small it would strongly suggest 
that the ss pairs in the proton are short-lived and play no essential role in understanding 
the proton. These are, indeed, important and interesting issues. It should be noted that 
there already is significant evidence that the non-strange sea is strongly interacting in the 



proton since the u/d radio is x -dependent [ p^ , p^ . However, it could well be that strange 
quark fluctuations are significantly shorter lived and interact more weakly in the proton 
than do the non-strange fluctuations. 

3. Heavy ion collisions and small-x physics 



3.1. Two main motivations 

There are at least two strong arguments for studing relativistic heavy ion colhsions. The 
first, of course, is that such collisions offer the possibility of producing, at least temporarily, 
a new state of matter, the deconfined quark-gluon plasma. The second motivation is to 
study high field strength QCD. At the very early stages of a heavy ion collision, well 
before equilibration occurs, very high values of the QCD field strength, F" , are reached. 
What the properties of such a system are, and how large F" can actually become, are 
fascinating questions which theorists have only recently began to investigate. 

3.2. Lessons from QED? 

In quantum electrodynamics electric fields, having a coherence length and lifetime 

— * 2 — * 

greater than 1/me, have a maximum attainable value, \E\ ~ ^^. Larger values of E 
are immediately shielded by the copious creation of e+e" pairs. Thus, crudely speaking, 
one can say that producing -^E"^ > 1/a will result in a breakdown of the QED vacuum. 
What about a hot QED plasma? How big are the field strengths in such a system? The 
question is relatively easy to answer. Suppose we look over a region of radius Vc and ask 
what is the size of the electric field coherent over a sphere that size. The answer is given 
by counting the number of photons in the plasma having \k\ < l/vc- Thus, 

rtE^ ex rljd'kQ{l/r, - \k\)—l—^ (22) 

Since [e^k/'^ ~ 1]~^ ~ "^ l^k ^^^ small \k\ it is clear that r^E"^ grows as r^ grows and that 
the growth stops when the plasma frequency up oc eT is reached. Thus for r^ ~ (eT)~^ 

vtEl oc 1/e (23) 

and this is much less than 1/a the maximum allowed value in the vacuum. The value 
of the maximum field found in a plasma is lower than that of the vacuum because the 
plasma has many electrons having momenta much greater than eT which are effective in 
shielding electric field fiuctuations whose value is greater than that given by (23). 

3.3. Crude picture of early stages of a heavy ion collision 

Likely at RHIC energies and certainly by LHC energies semihard gluon production will 
dominate the transverse energy freed in relativistic heavy ion collisions. At the moment 
there is a semiquantitative understanding of the early stages of a heavy ion collision 
coming from semihard gluon production p^-p8|. Imagine an ion-ion collision in the center 



of mass frame. Just before the collision the wavefunction of each of the ions has many 
small- X gluons. During the collision large numbers of these small-x gluons are freed by 
elastic gluon-gluon scattering as illustrated in Fig. 2. The freed gluons which dominate the 
produced transverse energy are expected to have p±_ ~ IGeV at RHIC and p± ~ 1 — 3GeV^ 
at LHC. At LHC the gluons are clearly within the hard scattering regime while at RHIC 
the hard scattering contribution should give a reasonable estimate of the freed transverse 
energy at early times after the collision. Rough estimates indicate that[P^-p7|] 



—^ ^ 1 TeVat RHIC (24) 

dy 



and 



dy 



15TeVat LHC 



(25) 



with the 1 TeV at RHIC coming from 10^ gluons and the 15 TeV at LHC coming from 
5x10^ gluonsp7| 




Figure 2. The dominant process in estimating gluon production in ion- ion colhsions 



3.4. A closer look at the early times after an ion-ion collision 



The McLerran-Venugopalan model |2^ is a nice framework within which to look more 
closely at the early stages of a high energy heavy ion collision. While the McLerran- 
Venugopalan model cannot be expected to be correct in the details of a heavy ion collision 
it should be a reasonable guide as to how the reaction proceeds. One begins by looking 
at the distribution of valence quarks in a high energy ion. The valence quarks are found 
in a Lorentz-contracted longitudinal disc of size Az = 2R ■ — where m is the nucleon 
mass and p the momentum per nucleon of the ion. For our purposes we consider Az = 
so that one can imagine the valence quarks having a two-dimensional number density of 
quarks per unit area 



ng{b) = Qp^R^ - 62, 



(26) 



where p is the normal nuclear number density, which we assume to be constant in the 
nucleus, while h is the impact parameter measured from the center of the nucleus in 
a direction perpendicular to the direction of motion of the nucleus. The situation is 
illustrated in Fig. 3. 

The valence quark density, given by (26), is the source for soft gluons corresponding to 
the Weizsacker- Williams field of Uq pO| , |3T[] . The color charge at a given impact parameter 
comes from a random addition of the color charges of each of the valence quark at that 
impact parameter. A single quark gives gives ^§^^n% gluons at scale Q^ per unit rapidity 
so that one expects the number density of gluons per unit area and per unity rapidity to 
be 



n,(6,Q)=n,(6)-^£nQ7/i^ 



TT 



(27) 
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Figure 3. The spatial distribution of valence quarks in a high energy heavy ion 



in the "additive" Weizsacker- Williams approximation. (If /i ^ lOOMeV (27) works rea- 
sonably well for a single proton.) However, a more careful application of the non-Abelian 
Weizsacker- Williams approximation shows that there can never be more than 1/a gluons 



occupying the same transverse area |30|j32|| . Thus 



n„ 



\b,Q) 



N^ 



1 



1 



naN, 7r(2/Q)2 






(28) 



where we take the area of a gluon at scale Q to be 47r/Q^ and we have inserted the color 
factors to agree with the results of Refs.30 and 32. Thus so long as (27) is less than (28) 
it should be a reasonable picture of the gluon distribution in a large nucleus. When Q^ 
is small enough so that (27) is greater than (28) one says that the gluon distribution has 



saturated|Q with n™"^ being the saturated distribution. We can make (27) a little less 
model dependent by identifying ^sc> inQ'^/^'^ with xG{x, Q"^) the gluon distribution in 
a nucleon. In that case (27) becomes, using (26), 



Ugib, Q) = 2pVW^¥ xG{x, Q2). 
Equating (28) and (29) gives 



Q] 



Nr 



'at 



in 



Ar2 

c 



-apVR^ - &2 xG{x, Q 



2 ) 

sat) 



(29) 



(30) 



as the saturation momentum. For Q"^ values below Ql^^ (28) should be the gluon distri- 
bution in the nucleus, while for Q"^ > Ql^^ (29) should be appropriate. The Q^^^ value 
given by (30) at 6 = is very close to that given long ago in Ref.24, differing only by a 
factor 3/4 and is the same as given in Refs. 30 and 32 though the discussion here has 
been much simplified. It is likely that the value of Qiat ^^ RHIC is about | — IGeV"^ or 
so while the value at LHC may well be in the Q^^^ ~ 2GeV'^ region at 6 = 0. 
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There is pretty good control of the ion's wavefunction in the McLerran-Venugopalan 
model, however, so far no one has succeeded in doing a realistic calculation of the "scat- 
tering" to determine which gluons are freed at early times. One would guess that all 
gluons having transverse momentum below Q^at would be freed while not so many of the 
gluons above Qsat become free. A good calculation of the gluon interactions to determine 
exactly which gluons are freed is a key calculation for progress in understanding the early 
stages of heavy ion collisions. 

3.5. What about the proton wavefunction at small x? 

Is the picture of the small-x gluon distribution in a nucleus unique to large nuclei or 
can the same picture apply to the proton? Saturation, or a maximum field strength, 
comes about when many gluons are available to occupy the same region of phase space 
in a wavefunction. In a large nucleus the source producing those gluons is the large 
number of valence quarks. In a proton a similar phenomena can occur at extremely 
small values of x, where now the large rapidity interval, y = in 1/x, available for gluon 
evolution can lead to high gluon number densities. Theoretical control over this very 
small-x regime is not perfect yet, and the importance of BFKL or x-evolution is not yet 
clear in the HERA regime, so one must turn to phenomenology to see if there is evidence 
for saturation effects. The answer is, perhaps, there is some evidence from low-Q^ very 
small-x structure function data at HERA. 




} 



quark-antiquarkpair 



Figure 4. Viewing the structure function of the proton in terms of the scattering of a high 
energy quark-antiquark pair. 



Let me briefly describe a useful way to view structure functions for the purposes of 
seeing if saturation effects are present. It is convenient to choose a collinear frame where 
the proton momentum, P, and the virtual photon momentum, q, take the form 
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Figure 5. The Caldwell plot of HERA data. 



with P/q very large and q/Q large but fixed independently of 1/x = ^^ Thus all x — 
evolution is included in the proton's wavefunction and the scattering can be viewed as 
illustrated in Fig. 4. Before the collision the virtual photon splits into a quark-antiquark 
pair which then interacts with the gluon field of the proton to produce an inelastic collision. 



A. Caldwell |3^ suggested looking at Q"^ — i^^a as a useful way of seeing very small-x 



eQ2 



effects at moderate Q^ . This is useful because 



where aqq is the cross section of the qq pair to interact with the proton. The transverse 
coordinate separation of the pair is proportional to IjQ. If gluons are reasonably dilute 
in the proton, so that the parton picture applies, one has 



a ^, 2 



o-qq ex — xG{x, Q ) (33) 
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Figure 6. Caldwell's compilation of points from fixed target experiments 



with G the gluon distribution in the proton. The 1/Q^ on the right-hand side of (33) 
reflects the diluteness of gluons in the proton. On the other hand, if gluons are packed 
densely in the proton one reaches the unitarity limit and cross sections become geometric 
in which case 



aqg oc tiRq 



(34) 



where Rq is the radius of these proton. Using (33) and (34) in (32) we see that a signal for 
saturation is that Q^^m should grow linearly in Q^ . HERA data is shown in Fig. 5 while 
fixed target data is shown in Fig. 6. Data are plotted as a function of x on the lower axis 
and as a function of Q^ on the upper axis. The turnover in the HERA data at Q^ ~ 2, 
and the lack of such a turnover in the larger x fixed target data may be the indication 
of parton saturation for all transverse momenta less than about 1 — 1 l/2GeV^[^5]-|37 . 



3.6. Why is equilibration nontrivial? 

We have seen that there are a great many gluons in the wavefunction of a high energy 
heavy ion and that many of these gluons are likely freed in a central ion-ion collision. 
With the large gluon densities present it would seem to be a simple matter to reach 
equilibration. However, the gluons are located in a single layer, perpendicular to the axis 
of collision. Thus immediately after the collision the gluon density will begin to rapidly 
decrease and it is not certain, a priori, whether equilibration will set in before the system 
falls apart. There are some pretty good indications from Monte Carlo calculations [^^-^ 



that kinetic equilibration does indeed occur at RHIC and LHC energies, however, it would 
be nice to see more clearly, in analytic calculation what the issues are and how certain it 
is that equilibration occurs. In what follows I give a very simplified version of a criterion 
for equilibration in the McLerran-Venugopalan model. 

We imagine a head-on collision of two heavy ions each of which has gluons in its wave- 
function saturated up to a scale Q. We suppose that during a time on the order of XjQ the 
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gluons having k± < Q, in the central unit of rapidity are freed while those having k± > Q 
are virtual fluctuations only and are not freed. Then follow a particular freed gluon in the 
central unit of rapidity. (We may suppose that kj_ of the gluon is k± ~ Q, since there is 
little phase space for k± << Q and freed gluons having k± » Q are few.) If this gluon 
has a collision with momentum transfer on the order of Q it has gone a significant way 
toward equilibration. Thus as an equilibration criterion we take 

dt , , 

where tq = l/Q, A(t) is the mean free path and r is some time which should obey r < 
R. (If the gluon does not experience a hard scattering in r < i? it is unlikely that a 
later scattering will occur because the expansion becomes 3-dimensional rather than 1- 
dimensional.) Now 

(36) 



\ /j.\ ~ Pgy'^J^gg^gg 




where 




Pg- 9it) = ^ = 


(iv2 - i)g2 


and 




"gg^gg v ^ / / a 


T2 T\n2 



(37) 



(38) 

In using (37) we neglect the b-dependence of the saturation momentum. In (38) we have 
integrated ^ = i^f (^Jl\y2 over t from -t = Q^ to oo. a in (37) and (38) should be 
evaluated at Q^. Thus (35) becomes 

1 = ^in t/to (39) 

TT 

or, using a(g2) = ^^^^ with b = ^^=^, 

r = i(^)2- b/N.-\ (40) 

Since Q^ cc R from (30) equilibration will occur if 

71 b 1 , ^ 

^-^<l (41) 

which is the case. It is reassuring to see that the equilibration seems to be met in the 
McLereran-Venugopalan model, but it is , perhaps, a little disturbing that it seems to be 
accidental and not a fundamental property of high energy heavy ion collisions. 
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4. Perturbative versus nonperturbative QCD 

Full nonperturbative calculations of physical quantities are desirable but hard to get 
in QCD. In hard probes of quark and gluon structure in hadrons the operator product 
expansion separates hard from soft physics where the hard part is calculable perturbatively 
and the soft part can be parametrized in a general way. The operator product expansion 
takes the form 

Je"^^d'x j{x)j{0) ^2^^^ Oii?i(a(Q)) + ^O^S^{a{Q)) + ■■■ (42) 

where Oi is the leading twist term and O2 and the next- to-leading twist term with Ei 
and E2 coefficient functions. An interesting subtlety arises because Ei has a divergent 
perturbation series 



E,{a) = ^^i„a" (43) 

n 

n large n 

where (52 is the first coefficient of the /?— function. The ambiguity in the sum (40) due to 
the divergent perturbation series is of size l/Q^ and means that the separation of the Oi 
and 



t02 terms in (39) is not without ambiguity ||4l[ . 

For example, for the F3 structure function in neutrino proton scattering one commonly 

writes 

xF,{x,Q^)=xF^^{x,Q^) + ^ (44) 

with F^'^ the leading twist contribution and the h/Q"^ the higher twist term. Recently 
an interesting observation has been made by Kataev]^ and his collaborators that h{x) 
depends on the level to which the perturbation theory analysis of F^'^ is done. If F^^ is 
evaluated in a leading order renormalization group formalism a large h{x) emerges from 
the fit. If a next-to- leading order formalism issued for F^'^ a somewhat smaller but still 
substantial contribution from h is needed. When F^'^ is evaluated in a next-to-next- 
to- leading formalism very little room is left for h. This is shown in Figs. 7- 9 taken from 
Ref.9. I think there is an important lesson here. It is not clear how to separate leading and 
higher twist terms and the separation will depend on the level to which the perturbative 
calculation is done. 
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Figure 7. The results of the LO extraction of the x-shape of the twist-4 contribution h(x) 
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Figure 8. The results of the NLO extraction of the x-shape of the twist-4 contribution 
h(x). For comparison, the IRR-model prediction of Ref. [22], obtained using the NLO 
MRS parametrization, is also depicted. 
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Figure 9. The results of the NNLO extraction of the x-shape of the twist-4 contribution 
h{x). 
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